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ABSTRACT 

We present a study of hydrodynamic drag forces in smoothed particle simulations. 
In particular, the deceleration of a resolution- limited cold clump of gas moving through 
a hot medium is examined. It is found that the drag for subsonic velocities exceeds 
that predicted by simple physical approximations. This is shown to be a result of 
the hydrodynamical method which encourages the accretion of particles from the hot 
medium onto a shell around the cold clump, effectively increasing the radius of the 
clump. For sonic and supersonic velocities, the drag is shown to be dependent on the 
effective cross section of the clump. The consequences for cosmological simulations are 
discussed. 

Subject headings: hydrodynamics - methods: numerical - galaxies: clusters: general - 
galaxies: kinematics and dynamics - large-scale structure of Universe 



1. INTRODUCTION 

Numerical simulation has proved to be a powerful tool for understanding the dynamics and 
formation history of cosmological structure (Davis, Efstathiou, Frenk, and White 1992; Katz 1992; 
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Cen, Miralda-Escude, Ostriker, and Rauch 1994; Evrard, Metzler, and Navarro 1996). When N- 
body techniques are combined with a hydrodynamical method to fohow the baryonic component, 
further advancements can be made as now no additional assumptions are required to relate the 
baryonic material to the underlying collisionless matter. Such a combination has proved to be both 
a popular and powerful technique (Hernquist and Katz 1989; Evrard 1990; Katz and Gunn 1991; 
Cen 1992; Navarro and White 1993; Steinmetz and Miiller 1993; Ryu, Ostriker, Kang, and Cen 1993; 
Bryan, Cen, Norman, Ostriker, and Stone 1994; Gnedin 1995; Steinmetz 1996). A common method 
of implementing hydrodynamics in N-body simulations is through the use of smoothed particle 
hydrodynamics (SPH) (Gingold and Monaghan 1977; Lucy 1977). The great strength of SPH is 
that it is a Lagrangian technique; individual particles act as tracers of the flow and better resolution 
can therefore be achieved in regions of high tracer density. This is a very desirable property for 
cosmological simulations where extremely large density contrasts are present and excellent spatial 
resolution needs to be achieved within relatively large volumes. 

In the standard model of cosmological structure formation, larger structures are formed via 
the amalgamation of smaller structures formed at earlier times in a hierarchical fashion. In any 
simulation modelling the formation of such structures in this hierarchical manner, there will be 
an abundance of objects formed at the resolution limit of the simulation. A particularly difficult 
problem is following these small, resolution-limited halos as they are subsumed into larger struc- 
tures as the hierarchy progresses. It has been found that much of the substructure in collisionless 
simulations is erased unless very high resolution is employed (Ghigna, Moore, Governato, Lake, 
Quinn, and Stadcl 1998; Moore, Lake, Quinn, and Stadcl 1999; Klypin, Gottliiobcr, and Kratsov 
1999). This is the classic overmerging problem (White 1976; van Kampen 1995; Summers, Davis, 
and Evrard 1995; Moore, Katz, Lake, Dressier, and Oemler 1996). 

With the inclusion of radiative cooling into the hydrodynamical model, the baryonic material 
can dissipate energy and become more tightly bound than the dark matter, allowing individual 
objects to survive the hierarchical process that builds up large structures (Summers et al. 1995; 
Prenk, Evrard, White, and Summers 1996). Klypin et al. (1999) argues that these objects are not 
likely representative of bound gaseous structures in nature, since they are lacking the dark matter 
halos observed to dominate similar structures (Persic, Salucci, and Stcl 1996). As well, the dynamics 
of this set of objects may now be affected by gas forces; as Frenk et al. (1996) demonstrates, the 
viscous drag due to the passage of a galaxy-size object through the hot intracluster medium can 
be sufficient to cause much of the material to spiral into the centre. 

An understanding of the numerical and physical processes involved as resolution- limited gaseous 
structures move around within the environment of larger halos is essential to understanding the 
results of numerical simulations. This includes the effect on the global properties of the larger 
halos. Velocity bias, where galaxies have a smaller velocity dispersion than that predicted by virial 
arguments due to the exchange of orbital energy to the collisionless material is sometimes observed 
in simulations (Evrard, Summers, and Davis 1994; Frenk et al. 1996) and sometimes not (Katz, 
Hernquist, and Weinberg 1992) (However, Summers et al. (1995) claims these results are strongly 
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afFected by the choice of galaxy tracer). Spatial bias, in which the galactic population is more 
or less centrally concentrated than the dark matter, again due to exchanging orbital energy with 
the collisionless material, is another effect sometimes observed in simulations of the formation of 
galaxy clusters (Cen and Ostriker 1993b). This can lead to incorrect estimates of the underlying 
density profile. Observationally, velocity and spatial bias have serious effects on mass estimates for 
clusters (Fabricant and Gorenstein 1983; Jones and Forman 1984; Girardi, Giuricin, Mardirossian, 
Mezzetti, and Boschin 1998) and hence on such fundamental quantities as the baryon density and 
the amount of matter in the Universe as a whole (White and Fabian 1995; White, Jones, and 
Forman 1997). Mass segregation, in which larger galaxies arc observed to lie, on average, closer 
to the centre of clusters than smaller galaxies, is another effect which occurs in many models of 
cluster formation (Frenk et al. 1996) . An understanding of the processes that affect both the rate 
by which the smaller structure is merged into larger halos and the final location of the deposited 
material are of fundamental importance. 

Any process which modifies the merger rate of small groups of galaxies in simulations also 
impacts the connection with observations and models of galaxy populations. It may have an 
important bearing on the physical processes that lead to the formation of giant cD galaxies. In 
addition, there has been much interest in the morphology-density relation in which elliptical galaxies 
are observed to be much more common in rich environments than spirals (Dressier 1980), a result 
supported by simulations (Elizondo, Yepes, Kates, and Klypin 1999; Cen and Ostriker 1993a). The 
merger of spirals to form ellipticals is currently a very popular model (White and Rees 1978; Cole 
1991; Lacey and Silk 1991; White and Frenk 1991) and the rate of this process will be determined 
by the detailed dynamics of galaxies within clusters. The merger rate of galaxies within halos is also 
an important parameter for the semi-analytic models of galaxy formation (Kauffmann, White, and 
Guiderdoni 1993; Cole, Aragon-Salamanca, Frenk, Navarro, and Zepf 1994; Baugh, Cole, Frenk, 
and Lacey 1998). These rates are often taken from low resolution SPH simulations similar to those 
we discuss here. 

The effect of drag on clumps of matter is significant in simulations of other astrophysical 
scenarios. SPH has been used to simulate the feeding of active galactic nuclei (Artymowicz and 
Lubow 1996; Sponholz and Molteni 1994), starbursting behaviour within massive galaxies (Heller 
and Shlosman 1994), formation of interstellar clouds through the merger of clumps (Kamaya 1997; 
Bliattal, Francis, Watkins, and Whitworth 1998), and, more recently, even planet formation (Nel- 
son, Benz, Adams, and Arnett 1998), all situations in which the merger rate of small objects is of 
vital importance. 

In Thacker, Tittley, Pearce, Couchman, and Thomas (1998), tests of a variety of implementa- 
tions of SPH are described. The implementations varied in the method of force symmetrisation, the 
form of the artificial viscosity, and the presence of a shear-correction term in the artificial viscosity. 
They were compared by probing the behaviour of the SPH implementations in a series of scenarios 
relevant to cosmological simulations. Included in this series is the drag experienced by a cold clump 
travelling through a hot medium. It was concluded that the drag on the cold clumps is excessive 
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when SPH is used, independent of the implementation method. The reason for this behaviour was 
not explored, nor was the behaviour extensively characterised. This paper is an extension of the 
drag tests work described in Thacker et al.. We present here the 

In this paper we investigate the numerical processes which control the strength of the viscous 
drag of a hot medium on a resolution-limited cold clump of gas within a smoothed particle hy- 
drodynamical (SPH) method. This is done through the use of numerical simulations of the simple 
scenario in which a cold gaseous halo is passed through a tenuous hot, uniform medium. Three 
velocity regimes are explored: subsonic, sonic, and supersonic. Analysis is assisted by a compari- 
son of the actual deceleration observed for our clumps with that predicted from simple theoretical 
arguments. 

Our model clumps are intentionally simple. We study what are essentially resolution dominated 
spheres of material moving within a uniform density hot medium. Such objects arc typical of the 
small galaxies or groups of galaxies formed in cosmological simulations of structure formation where 
most of the objects have masses and physical extents near to the resolution limit of the model 
(Pearce, Jenkins, Prenk, Colberg, White, Thomas, Couchman, Peacock, and Efstathiou 1999). It 
is important to note that regardless of the resolution, these resolution-limited structures will form 
and play a role in the development of the larger structures. 

The paper is laid out as follows. In Sec. 2 we include a brief overview of the relevant SPH 
equations highlighting those factors which affect the drag. The model clump-halo systems are 
described in Sec. 4. The predicted deceleration for a galaxy which sweeps up all the material it 
encounters is derived in Sec. 5. A comparison between this and the actual deceleration found in 
our models for a variety of Mach numbers is presented in Sec. 6 and discussed in Sec. 7. 



Under the SPH formalism any scalar quantity at any position can be calculated by averaging 

the properties of the surrounding particles. This average is generated by weighting the contribution 
of each particle as a function of the distance from the designated point using a weighting function 
(or kernel) that is smoothly varying. 

The average of any scalar quantity, ^(r ), is given as; 



2. 



SMOOTHED PARTICLE HYDRODYNAMICS 



2.1. Force symmetry 




(1) 



where W{r ,h) is the kernel function and h is the smoothing length. The kernel function is generally, 
but not necessarily, spherically symmetric, has compact support and is normalised such that its 
spatial integral is unity. 
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The gradient of any scalar field can be calculated by using the gradient of the kernel function. 
Noting that the integral becomes a summation over N neighbouring particles for a finite distribution 
of particles, the gradient estimate is, 

AT . , , 

where rij is the local particle number density. 

In the above description, the smoothing length, h, has been taken to be a constant. However, 
to make the best use of the Lagrangian nature of the technique, this would be inappropriate since it 
would establish a minimum resolution scale far larger than the resolution permitted by the increased 
number density of particles in the dense regions. To allow for this, the smoothing length should 
decrease as the local number density increases to maintain an appropriate spatial resolution. For h 
to scale appropriately with the resolution it is normally taken to be a distance which encompasses 
some constant number of neighbouring particles, Nsph- 

Allowing h to vary leads to the question of which particle's smoothing length should be used 
in the SPH summations. There are two interpretations of SPH; the 'gather' and the 'scatter'. In 
the 'gather' interpretation, a local estimate of a property about a particle, is made by sampling 
out to some fixed distance and weighting each particle found by its separation from the centre, 
hence the smoothing length in the summation is hi. In the 'scatter' interpretation a local esti- 
mate is made using the sum of the contributions to that position from all the particles for which 
their smoothing length extends to encompass the position. Here, hj is used. In practice, neither 
interpretation is accepted but rather the appropriate smoothing length is some function of the 
individual ^'s, hij = f{hi,hj). The form of hij can be the arithmetic, harmonic or geometric mean 
of the ^'s, or some other form. If it is symmetrical, then f{hi,hj) = f{hj,hi). It need not be 
symmetric in the general case (for example, to calculate densities), but if it is to be used in the 
calculation of the force, symmetrisation will preserve force symmetry and conserve momentum. 
A further, equally valid, option is to average the kernels themselves; W(ri — rj,h) is replaced 
with {W{ri — r j,hi) + W{ri — r j, hj)) /2. This also preserves the force symmetry and conserves 
momentum. 

To demonstrate the utility of these symmetrisation methods, let us examine the basic hydro- 
dynamical force equation. The hydrodynamical component of the force, F, on each particle comes 
from the pressure gradient, VP, via 

VP 

F = -m—, (3) 

where m and p are the gas particle mass and density, respectively. To preserve force symmetry it 
is necessary to find an SPH approximation of the hydrodynamical force contribution such that Fij 
= —Fji where F^j is the force felt by particle i due to particle j. This is frequently achieved by 
using the identity; 

VP P P 

— = V- + ^Vp. (4) 
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Combining this with Eq. 3 we get, 




This is not yet symmetric. When the force on particle j is computed, the contribution from 

particle i, Fji, will not be the negative of the contribution Fij to Fi since VW{ri — rj,hj) ^ 
— Vty(r j — r j, /ij) if hi ^ hj. However, if hi and hj arc replaced with a symmetric combina- 
tion, hij, it is clear that the gradient will also be symmetric. For a spatially varying smoothing 
length, Monaghan (1992) suggests taking the gradient of some average of the kernels calculated 
symmetrically; 

Fi = -rmmj ^ ( 4 + 4 J ^^ir (6) 

j \Pj Pi J 

The 'average' kernel, Wjj, can use the kernel averaging scheme seen above or some average of the 
smoothing lengths for hij. That is, either 

VWij = [VW{vi-rj,hj) + VW{ri-Vj,hi)\/2 (7) 

or VWij = VW{ri — rj,hij) where hij is an average of the smoothing lengths seen before. The 
former approach is used in Steinmetz and Miiller (1993). 

The effective kernel shapes produced by two smoothing radii that differ by a factor of two are 

illustrated in Fig. 1. The effective kernel shapes differ most appreciably in the wing of the function, 
particularly between the /i-averaging and the kernel averaging schemes. For the massive cold clump, 
the increase in the wing for the kernel averaging scheme increases the effective cross-section of the 
clump and, consequently, increases the drag (see Sec. 6.1). 



2.2. Artificial viscosity 

Early on, it was realised that for SPH a viscosity force term would be necessary to damp the 
flow of gas (see, for example, Lucy (1977)). This is entirely to be expected; the gas is represented 
in the simulation by particles of macroscopic size, while viscosity is due to turbulence on all scales, 
including microscopic. The cascading of energy from large-scale motions to small-scale turbulence 
and eventually, via atomic interactions, into heat must be modelled by a viscosity term in the 
hydrodynamic equations. The adopted forms of the viscosity terms are essentially ad hoc. They 
are all designed to cover bulk viscosity and shock front dissipation (Monaghan 1992). In some 
cases, it is important to modify these forms to permit shear (Balsara 1995). All, however, convert 
convergent kinetic energy into thermal energy. 

It will be shown in Sec. 6.2 that the triggering of the artificial viscosity does affect the drag felt 
by the clump in the subsonic case. However, the form of the artificial viscosity is not significant to 
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the general behaviour of the drag, as is shown in Thacker et al. (1998). Interestingly, they report 
that artificial viscosity actually decreases the amount of drag experienced by cold clumps, since 
it reduces the penetration of hot particles into the cold clump. Without the artificial viscosity, 
the convergent flows experienced ahead of the clump are not damped, allowing particles to pass 
through the bow wave and approach closer to the cold clump. This increases the instantaneous 
pressure force between the hot particle and the cold clump leading to a collision more elastic in 
nature. Elastic collisions take away up to twice as much momentum as inelastic. With the artificial 
viscosity term, the hot particles are diverted into the flow around the clump, heading off any close 
approaches. 

3. NUMERICAL METHOD 

The implementation of SPH used here was kindly made available by Dr. R. Thacker and 
is described in Thacker et al. (1998). It is based on the N-body hydrodynamic code, hydra 
(Couchman, Thomas, and Pearce 1995). Further details about the code are available there. A 
more complete description of the methods for generating the symmetrical hydrodynamical forces 
arc contained in Thacker et al. (1998) as well as the implementation of the form viscosity. We use a 
Monaghan-type artificial viscosity. Symmetrisation of the force equation is provided by arithmetic 
averaging of the smoothing lengths unless otherwise noted. 

Since it was demonstrated in Thacker et al. (1998) that the general behaviour of the drag is 
not affected by the presence of a Balsara term (Balsara 1995), which reduces the viscosity in the 
presence of shear, one was not included in the implementation of SPH used. 

For the calculation of the gravitational forces, HYDRA uses an adaptive particle-particle, 
particle-mesh (AP^M) method. The minimum resolution is set by the gravitational softening 
length which is 75 kpc . 

All forces are periodic across the boundaries of the simulation cube. 

4. MODEL SYSTEMS 

The tests in this paper are all based on a model containing a cold clump moving through a hot 
medium. This clump may be allowed to slow or be forced to maintain a constant velocity, depending 
on the needs of the particular test. Our models attempt to mimic the conditions encountered as a 
resolution-limited cold clump passes through the halo of a larger structure. The models consist of 
a knot of 100 cold particles moving through a uniform distribution of 13000 hotter particles. The 
volume is 2 Mpc on a side and has periodic boundary conditions. 

The clumps of cold gas are essentially structureless spheres, lying at the resolution limit of the 
code since they are to represent the first (and most common) structures to form in a cosmological 
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simulation of hierarchical clustering which is the dominant source of material during the creation of 

larger structures (Kay, Pcarcc, Jenkins, Frenk, White, Thomas, and Couchman 1999). The gas in 
the clumps is cold (10^ K ) since the dense gas cools rapidly to this temperature through a radiative 
process. 

Within a typical cosmological simulation the halos of galaxy clusters are hot (< 10^ K) and 
more diffuse than the cold clumps because the hot halo gas is sitting close to virial equilibrium 
within the cluster potential well. This hot halo is observed through the X-ray emission (see Horner, 
Mushotzky, and Scharf (1999), for examples). Consequently, the hot media in our models are at 
temperatures of 10^ or 10^ K and are 1/100 as dense as the cold clumps. 

To cover a variety of infall speeds, we examine the deceleration of a knot of cold gas in three 
velocity regimes: Mach 2, Mach 1, and Mach Vs- The Mach 2 and Mach 1 tests differ in terms of 
the speed of the cold knot; the temperature of the hot gas is kept fixed. The Mach l^ test uses the 
same clump velocity as the Mach 1 test, but is performed in hotter gas. Table 1 gives the details 
of the cold clump and hot gas phases. 

The hot gas is in a glass configuration which was prepared from an initially random placement 
of particles and subsequently allowed to relax to a stable state. The cold clump was created by 
randomly placing stationary particles within a sphere of radius equal to the gravitational softening 
length then allowed to relax in isolation. The two systems were then combined. 

The Jeans length, Rj, for the hot gas phases is sufficiently large to ensure stability even in 
the presence of the perturbation from the cold clump. Consequently, dynamical friction should 
not be important. This conclusion was confirmed by passing a collisionless clump through the hot 
medium - it experienced negligible deceleration. Ostriker (1999) discusses dynamical friction in 
systems with hydrodynamics. 

The box length, 2 Mpc , was chosen so that the cold clump was well separated from its images 
(arising from the periodic boundary conditions employed) and would move across the box only once 
without encountering its own wake. 



The expected maximal rate of deceleration can be approximated by considering a disc sweeping 
through a medium, collecting all the material it encounters. This would represent the maximum 
expected rate of deceleration if dynamical friction and elastic collisions are ignored. Prom the 
equations of conservation of mass and momentum, the solution for the velocity, V, of such a system 
is found to be 



5. 



EXPECTED DECELERATION 




(8) 
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where ti/^ is a characteristic time scale corresponding to the time at which a mass of the hot gas 
equal to one half that of the disc has been swept up. It is given by 

M 

Here, M is the mass of the disc at the start, R is the radius of the disc, and pg is the density of 
the gas through which the disc is travelling. The disc starts with velocity at time t = 0. The 
expected acceleration is just the derivative of this expression, and at i = 0, 

dt ° M ^ ' 

Ruderman and Spiegel (1971) derives a similar expression for the contribution from bow shock 
drag to the deceleration of a galaxy passing through an intraclustcr medium. Its form differs in 
that it contains a factor for the coefficient of drag, Cd, which is approximately unity. 

For our models the values of M, Vq and pg arc given. For any decelerating system, the drag 
may be associated with Eq. 10, allowing the calculation of an effective radius, Reff- An estimate 
of Reff may be made as follows. For the clump particles, the smoothing length, hdump-, is set to 
the minimum allowable value, as determined by the gravitational softening length. This leads us 
to choose Reff = "^hciump- However, because of the symmetrisation inherent in the SPH method, 
the smoothing length of the hot particles is also important and it is some combination of this 
length with the smoothing length of the cold particles that is required. We will assume that the 
appropriate radius is that radius at which the clump particles contribute an equivalent amount to 
the SPH density calculation for a hot gas particle as the hot particle itself. That is, 

NclumpWiReff,hij) = WiO, hhot)- (11) 

Since Ndump » W{0, hhot), this implies W{Reff, hij) ~ which gives -Re// — 



6. RESULTS 

Using the model clump-hot-gas systems, a series of tests were performed to characterise the 
behaviour of the drag in three velocity regimes: subsonic, sonic, and supersonic. A description of 
the individual tests as well as their results are presented in the following subsections. The results 
of these tests are broken into those that pertain to the clumps in all of the velocity regimes and 
those that are peculiar to the subsonic regime. 



6.1. Viscous drag in three velocity regimes 



The behaviour of the drag that is common to the velocities spanning subsonic to supersonic 
is described here. First, we will demonstrate that the viscous drag for subsonic clumps is much 
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greater than that expected from theoretical arguments. For all velocity regimes, it will then be 
shown that the effective radius of the clump, as determined by the drag it encounters, scales with 
the averaged smoothing length, hij. It will be verified that as a consequence, the drag varies with 
the method of force symmetrisation. Since hij is also a function of number density contrast, the 
increase in drag in regions of higher density will be demonstrated to be partially offset by the 
decrease in the averaged smoothing length. The distribution of the energy deposition will also be 
examined. 



6.1.1. The measured drag versus the fiducial drag 

The fiducial drag, given by the analytic expression Eq. 10 amounts to a maximum amount of 
expected deceleration, presuming inelastic collisions and a non-rigid medium. A measured drag 
which exceeds the corresponding fiducial value for the clump would necessarily by excessive. Mea- 
surements of the drag on the clump as it travels at a variety of velocities compared with the fiducial 
amount delineate those velocities for which the drag is excessive. To this end, a series of tests was 
conducted which allowed the drag to be calculated for such clumps. 

The tests involved ploughing a clump through a uniform medium at various fixed velocities. 
The mean drag on the clump was measured from the force required to keep the clump moving at 
a constant speed. The drag, as it varies with Mach number, is illustrated in Fig. 2. The results 
indicate a linear trend of the drag with velocity for speeds greater than about Mach 1/2- 

If the fiducial force of drag given by Eq. 10 is assumed and the radius of the clump, R, is 
taken to be twice the combined smoothing length, hij, as detailed in Sec. 5, then any points that lie 
above their respective curves in Fig. 2 indicate Mach numbers for which the drag on the clump is 
excessive. This is clearly the case for velocities below Mach 1. The Mach number of the transition 
to the domain of excessive drag is similar regardless of the sound speed for these tests, despite the 
fact that the speed of sound varies by a factor of three. 

The results may be illustrated by another method. The measured drag force equated with the 
left hand side of Eq. 10 gives an effective radius, Reff- This radius, and consequently the effective 
cross section of the clump, is much larger than the estimated value of 2hij at slow speeds, which 
is shown in Fig. 3. Indeed, it is larger than 2hhot, the max;imum sphere of direct influence possible 
through SPH force calculations. 

At supersonic speeds there is an apparent convergence of the effective radius to a value less 
than the estimated 2hij. Indeed, it seems to be converging to hij, which is not surprising since 
the density of the bow shock will increase as the clump velocity rises while Eq. 11 assumes the 
hot particles are relatively isolated compared with the density of the clump. However, the relation 
between the effective radius and the Mach number, M. , 

^ = 1.66 ± 0.08 M -0-59±o.o4^ (12) 
hij 
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illustrated in Fig. 3 also fits well and it converges to zero for large M. . Thus, there is not enough 
information to argue for a convergence to any particular value. 

The effective radius of the clump may be made consistent with the estimated value of 2hij if 
we include a coefficient of drag in Eq. 10 with a value Co = 0.5. 

It is notable that the effective radius is not a function of the speed itself, but of the speed 
scaled by the local speed of sound. This emphasizes the velocity scale invariance of the results 
presented in this paper. 

The excessive drag at subsonic speeds manifests itself in Fig. 2 by the excessive drag over the 
fiducial force of drag and in Fig. 3 as the rapidly increasing effective radius towards lower speeds. 
We will return to this point in Sec. 6.2 where the cause of this phenomenon will be explored. 

6.1.2. The scaling of the drag with the averaged smoothing length 

If the effective radius of the clump is truly set by the radius of the combined smoothing length, 
hij, as is argued in Sec. 5, then the drag on the clump should scale with this length. To test this the 
clump was passed through three media of the same physical density, but differing number density, 
n/iot, of hot gas particles. Since the smoothing length, hhot, scales with Ujl^^, this varies the mean 
smoothing length of the hot gas and consequently the combined smoothing length, hij. The test 
was done for number density increments of two and five times. In these tests, which were performed 
for the velocity regimes of Mach Va; 1; and 2, the chimps were permitted to slow instead of being 
held at constant velocity as was done for the tests in the previous section. 

For each of the tests, an effective radius, Reff, was determined by fitting the velocity of the 
clumps to the form of Eq. 8. The results, given in Table 2, support the hypothesis that the 
effective radius does indeed scale with the averaged smoothing length, hij. In each velocity regime, 
the effective radius is approximately a constant multiple of the averaged smoothing length (fifth 
column). The correlation is not strong enough to be distinguished from scaling with just the 
smoothing length of the hot gas, hfiot (fourth column). It is clear, however, that the effective radius 
does not scale with the smoothing length of the clump particles (third column). 

Though the effective radius is a similar factor of the averaged smoothing length within each 
regime, the factor is different between each regime; Reff — 3 for the Mach l^ tests, 2 and 1.5 for 
the Mach 1 and Mach 2 regimes, respectively. This trend is the same as outlined in Sec. 6.1.1 and, 
in particular. Fig. 3. These scaling coefficients may be compared with the estimated factor of 2 
made in Sec. 5. Only the deceleration of the Mach 1 clump matches, while the effective radius of 
the supersonic clump is slightly less than the estimate. The result for the Mach Vs clumps further 
emphasizes the excessive amount of drag found in the subsonic regime. 
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6.1.3. The scaling of the effective cross section with local density 

The results of the drag tests involving the deceleration of a clump through media of varying 
number density (Sec. 6.1.2) indicate that the drag felt on a clump travelling in the outer halo of a 
gravitationally bound system will be increased due to the low number density of particles in this 
part of the halo. As the clump travels into denser parts of the halo, the apparent cross section will 
be reduced due solely to the increased number density of the surrounding particles. This should 
offset to a degree the increase in drag that the clump would be expected to feel due to the increase 
in the physical density. The force of drag in Eq. 10 is related to the effective radius and the local gas 
density via F^rag —R^ffPhot- Since we have the relations, Reff oc hij oc n^^J'^ , then Reff oc Pf^^j.^ 
may be derived. This implies that the drag force for the clump scales with the local gas density as 

Fdrag OC —Phof 

This result assumes hij oc /i/jot, an approximation that is always true for force symmetrisation 
by kernel averaging and true for symmetrisation by arithmetic averaging only for hhot » hdump- 
For harmonic averaging of the smoothing lengths, this reduction in drag as density increases is only 
partially effective for small hhot and disappears for hhot >> hdump- Geometric averaging leads to 

2 / 

the different relationship of F^rag oc —ph^v which is closer to the form normally expected. 

To verify that this offsetting phenomenon occurs, a set of tests similar to those in Sec. 6.1.2 
were performed but with the mass resolution held constant while the number density is varied. 
Hence, the clump was ploughed through material of both different physical densities as well as 
number density. Since the velocity of sound varies only with the temperature of the gas, varying 
the mass density does not change the Mach number. 

The measured deceleration for the clump in each test was then fit to the analytic form given 
by Eq. 8 in the same fashion as described in Sec. 6.1.2. This analysis produced an effective radius 
which could then be compared with the smoothing radii of the clump and the hot particles as 
well as the averaged value. If the effective radius still scaled with the averaged smoothing radius, 
/ijj, then the offsetting phenomenon would be verified. The tests and analysis were done for hot 
media with physical densities of 1, 2, and 5 times that employed in Sec. 6.1.2 for the velocities 
corresponding to Mach Vs, Mach 1, and Mach 2. 

The results of the analysis, summarised in Table 3, confirm that as the density contrast between 
the clump and the background density decreases, the effective radius still scales with /ijj and, as a 
consequence, the force of drag scales with pj^^^. In Table 3, the scaling of the cross-section with the 
averaged smoothing radius, /i^j, (fifth column) is very clear and distinguishable from scaling with 
either of the other two smoothing radii, hdump or hhot (third and fourth columns). This is true 
independent of the velocity regime. 
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6.1.4- The method of force symmetrisation 

Since the effective radius of the clump scales with the averaged smoothing length, hij , then it 
is to be expected that the method of force symmetrisation, which requires the use of an averaged 
smoothing length, should affect the drag. Verification of this is significant, since the method of 
force symmetrisation is something that is chosen during the coding of SPH. 

For the method of force symmetrisation provided by using a common value, hij , of the smooth- 
ing lengths calculated by arithmetic, harmonic, or geometric averaging, it is easy to show that 
arithmetic averaging will always produce the largest value of hij while harmonic averaging leads 
to the smallest. Averaging the kernels themselves does not lend itself to such a comparison, since 
the effective value of hij would be dependent on r. However, using the argument to estimate 
Reff ~ 2hij outlined in Sec. 5, for kernel averaging we have hij ~ hhot, which is a maximum. 
Figure 1 verifies that at large radii, averaging the kernels indeed leads to the larger weighting and, 
hence, degree of interaction. 

For those model systems in which the ratio of hhot to hdump is 2.7 times, the differences between 
the method of force symmetrisation should be quite appreciable; averaging the kernel produces an 
estimated radius of interaction almost twice that of harmonic averaging of the smoothing radii. 
This is consistent with the results of Thacker et al. (1998) in trend, if not magnitude, in which it is 
noted that averaging the kernels leads to the largest amounts of drag at velocities less than Mach 1 
while harmonic averaging of the smoothing lengths leads to the least. 

As a test, the drag on a constant- velocity Mach 1 clump was measured while varying the method 
of symmetrisation. The implementations used kernel, arithmetic and harmonic averaging schemes. 
Drags of 1.4 x 10^^, 1.0 x 10^^, and 0.6 x 10^^ dyne were measured respectively. The test was not 
performed for geometric averaging of the smoothing lengths since this method produces a value for 
hij intermediate between arithmetic and harmonic averaging. These results confirm the trend that 
larger values for hij produce a greater drag, with the kernel averaging scheme producing the largest 
mean drag, 2.4 times larger than that produced with harmonic averaging of the smoothing lengths 
and 1.4 times larger than that found when arithmetic averaging is implemented. However, since 
the drag should scale with R^ff-, kernel averaging should produce four times as much drag, instead 
of the 2.4 measured. This indicates that there are secondary effects produced by the method of 
symmetrisation that, for these tests, offset the direct connection between hij and drag. However, it 
is clear that the choice of method for force symmetrisation significantly affects the amount of drag 
incurred. 



6.1.5. Energy deposition 

The energy lost by the deceleration of the clump goes primarily into heating the surrounding 
gas. Since this gas is already hot, this extra thermal energy is absorbed without any significant 
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change in the temperature of the gas. The total energy of the clump is less than 1% of the total 
thermal energy of the hot gas. The spatial distribution of this deposition of energy into the hot gas 
can be examined by looking at the change in temperature of the gas particles. This is illustrated 
in Fig. 4. As expected, a bow shock is formed in the sonic and supersonic cases. On account of the 
large filling factor of the shock, it contains a large fraction of the deposited energy. Energy is also 
deposited in the wake of the clump. 

For the subsonic clump, a transitory pulse is created at the start of the simulation which 
quickly dissipates. Subsequently, most of the energy is deposited locally and carried along with the 
clump. 

6.2. Drag in the subsonic regime 

It was established in Sec. 6.1.1 that subsonic clumps have inappropriately large effective cross 
sections. This is demonstrated by Tables 2 and 3 and Figures 3 and 2. Whereas the sonic and 
supersonic clumps always decelerate less than the fiducial rate, the subsonic clumps experience a 
greater drag than expected. We show in this section that the clump accretes a layer of hot gas 
particles into a shell at the effective smoothing radius, 2/iij, and it is this shell that enhances the 
drag. Furthermore, it is demonstrated that this is a by-product of the SPH pressure calculation 
in the vicinity of a compact collection of cold particles coupled to the artificial viscosity which 
damps rapid accelerations. In this configuration, the calculated force on a hot gas particle due to 
gas pressure causes oscillations in the hot particle's velocity which triggers the viscosity term of 
the hydrodynamic forces, bringing the hot particle to rest relative to the moving clump. In the 
subsonic regime we expect the clump to accrete some material as it moves through the hot medium, 
but it is the increase in cross section due to this shell of accreted particles that causes the clump 
to slow excessively. 

6.2.1. The accretion of remoras 

Using the data from Sec. 6.1.2 in which a cold clump was allowed to decelerate in a uniformly 
dense medium of hot gas, the behaviour of the drag on a subsonic clump may be characterised. 
The subsonic clump initially follows approximately the same deceleration curve as given by Eq. 8. 
The transition to a more rapid deceleration occurs as particles are accreted from the hot gas phase. 
Accretion of these particles is evident in Fig. 5 which displays the flow of particles about the 
cold clump in the reference frame of the cold clump at various Mach numbers. No hot particles 
are decoupled from the flow of hot particles in the Mach 1 and Mach 2 simulations while 44 are 
decoupled in the Mach Vs run. 

The particles picked up from the hot medium will be referred to as 'remoras'. Since remoras 
are also found about the subsonic clumps which are not permitted to decelerate (as in Sec. 6.1.1), it 
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is clear that the accretion does not occur only as the clump slows to almost stationary, but occurs 
while the clump still has a significant velocity. 

The remoras form a narrow shell about the clump at a radius equal to the effective smoothing 
radius, 2hij. Note that in this section we are using simple arithmetic averaging for the combined 
smoothing length; that is hij = (hhot + hdump)/'^- This shell is indicated by Fig. 5 and is verified 
in Fig. 6, which plots the 2-point correlation function about the clump particles. The correlation 
function of the hot gas particles (dashed line in Fig. 6) shows that the particles have preferential 
spacings typical of a glass distribution, but not at the distance of 2hij. This association with the 
radius hij is strong evidence that the SPH calculations are at the root of this phenomenon. 

The tight span of radii in which the remoras are permitted to come to rest is demonstrated 
by the narrowness of the peak in the correlation function at hij, witnessed in Fig. 6. The jump in 
the correlation function is two orders of magnitude and is located just beyond 2hij with a width of 
15kpc or about 10% of 2/ijj. 



The presence of this shell of remoras about the cold clump is clearly a by-product of the 
hydrodynamical calculations; it is not seen in collisionless simulations and the radius of the shell is 
coincident with the averaged smoothing radius. Here, we will examine the behaviour of the SPH 
force calculation for a hot gas particle in the vicinity of a dense knot of cold gas. It will be shown 
to lead to a process that traps the hot particles at the radius, 2hij, seen for the remoras in the 
simulations. 

In SPH, the pressure forces on an individual particle, Fj, are summed symmetrically over all 
neighbouring particle pairs in the manner of Eq. 6 which can be rewritten 



This utilises SPH's ability to estimate gradients of local scalar quantities by scaling the local values 
of these quantities (the pressure, P ~ pT, here) by the gradient of the kernel function, W. This is 
normally very effective when the local gradients in the temperature and, more significantly, the gas 
density, p, are small. While the temperature is a parameter accumulated over time, p is a quantity 
re-estimated each timestep. Hence, in the presence of the steep density gradients surrounding cold 
clumps, this value can fluctuate significantly between iterations. If a hot gas particle is within 
twice its smoothing length, 2hhot-, of the centrally-concentrated cold clump, then the estimated gas 
density will be significantly larger (1.5 x, here) for the particle than otherwise. 

It is informative to determine what the summation of the pressure forces is for a hot gas 
particle between 2/ijj and 2hhot from a cold clump. There will be two regimes: one in which the 
gas particle is at the far end of this distance, and the other when it is near 2hij. In the first case. 



6.2.2. SPH forces for a particle near 2hij 




(13) 
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the hot gas particle will not be surrounded by an isotropic distribution of particles. There will be 
a void on the side of the particle facing the cold clump. The 2-point correlation function of the 
particles in the vicinity of the cold clump (Fig. 6) demonstrates the existence of this void which 
occurs because particles are repelled by pressure forces close to the cold clump. The particles in 
the cold clump are centrally condensed within a radius much smaller than their SPH smoothing 
radius which is, in turn, half the gravitational softening radius (the minimum value we allow it to 
take). Surrounding this clump is a void with a radius slightly smaller than the smoothing radius 
used in the SPH calculations when dealing with a pair of hot and cold particles. Hence, the average 
gradient will force the particle toward the cold clump. 

Once the particle is within 2hij of the cold clump it suddenly sees many more than Nsph 
neighbours so its smoothing length drops while its density estimate jumps. For a very diffuse hot 
medium essentially all the neighbours of the particle are now in the cold clump and the particle is 
repelled strongly by pressure forces. The kick the particle receives from the interaction with the 
cold clump is short lived, and soon the unbalanced pressure force from the surrounding hot gas 
forces the particle back towards the clump starting the cycle again. 

A hot gas particle near a cold clump thus oscillates (or perhaps, more correctly, bounces) just 
beyond the radius, 2hij. This oscillation is a consequence of the form of the SPH force calculations 
coupled with its treatment of neighbours and, as such, will be a characteristic of any implementation 
of SPH. 

The path of a particle that is just grazed by the cold clump as it passes is shown in Fig. 7. The 
path, in the reference frame of the cold clump, shows the scries of accelerations described. What is 
even clearer is how these accelerations conspire to trap the particle at a radius approximately half 
way between the smoothing radii for the clump and the hot gas. 

In Fig. 8, which plots the hydrodynamic forces for the hot gas particle in a direction perpen- 
dicular to the direction of motion of the clump, the alternating directions of the forces are clearly 
illustrated. In this figure, the density of the hot particle, as estimated by the SPH density calcu- 
lation, is shown at the bottom. It clearly delineates the two regimes described previously; when 
the density is low, the particle does not 'see' the clump whereas when it is high, it is under the 
influence of the clump and feels a strong force away from it. 

The oscillations between these two regimes arc important as they trigger the artificial viscosity 
which drags the particle into moving with the clump. As the viscosity is only triggered if two 
particles are converging, the particle is interacting viscously alternately with the clump and the 
hot gas. The viscosity attempts to damp the oscillations, which can also be seen in Fig. 8. The 
attempt to damp the oscillations has the by-product of turning the particle into a remora; as Fig. 9 
shows, though the pressure is still the dominant force, the accumulated momentum in the direction 
of travel is controlled by the contributions from the viscosity term. 



6. 2. 3. The significance of a compact core 



The previous description of the problem indicates that the compact nature of the clump exac- 
erbates the problem of the creation of a shell of remoras; the hot particles are influenced by either 
all or none of the clump particles causing large swings in the pressure forces. The effect may be 
mitigated if the hot particles come under the inflTicnce of the cold particles gradually, instead of 
abruptly as they approach the clump. This may be tested by having the clump particles occupy a 
larger volume within their smoothing radius. Presented here is the method by which this was done 
and the resulting effect on the drag. 

The cold clump particles may be forced to form a more extended clump by modifying the 
form of the gravitational potential to produce a repulsive force at close distances. The unmodified 
gravitational force expression is a spline approximation of a Plummer-softened 1/r potential. This 
was replaced by a gravitational potential of the form 

(14) 

^ x+1 ^ ' 

where x = r/tg with being the gravitational softening parameter. Experimentation found that 
the parameters a = 0.2 and n = 2 suitably dispersed the clump, increasing the mean interparticle 
spacing by a factor of 100, but keeping the clump particles within twice the gravitational softening 
length. 

The influence of the antigravity factor extends to beyond twice the gravitational softening 
length to compensate for the strong compressional forces exerted by the pressure of the surrounding 
hot gas. As such, this form does not permit a stable clump to exist in isolation. This also does not 
permit a simple drag experiment; by the time the clump has expanded to a new stable configuration 
within the hot gas, the clump has already decelerated appreciably. Instead we again hold the 
velocity of the cold clump fixed as it progresses through the hot medium. This allows a stable 
configuration to develop and the drag on the clump to be measured. 

The drag for the compact clump was measured to be 6.0 x 10^^ dyne while the drag on the 
more diffuse clump was 4.2 x 10^^ dyne . The diffuse clump experiences almost Vs less drag. 

This factor is remarkable, considering that the diffuse clump has a larger cross-section. The 
cross-sections may be measured by convolving the clump particle positions with their effective 
smoothing radius for interaction with the hot particles, /ijj, to produce a column density distribu- 
tion. Since the hot gas particles are influenced strongly at even small column densities, we take 
10% of the maximum in the column density distribution for the compact core to calculate the total 
effective cross section of the clumps. The extended clump is found to have a cross-section that is 
60% larger than that of the compact core. 

Taken together, that the diffuse clump experiences Vs less drag while having a cross section 60% 
larger implies that the effect of the abrupt transition between the cold clump and the hot particles 
caused by having a very small central clump doubles the drag. This is a clear demonstration of the 
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significance of a compact core in enhancing the drag. 

Thacker et al. (1998) notes in a similar set of tests a difference in the clump sizes depending 
on the viscosity form adapted. It reports that a Monaghan-type artificial viscosity, as used in these 
tests, produces a clump that is two to three times more extended that those that use a V • v 
form. It should be emphasized that the mean intcrparticle spacing of the cold clump particles was 
required to be expanded by a factor of one hundred before the drag was diminished. Hence, the 
increase in drag is not something that can be easily fixed by modifying the implementation of SPH. 
This is particularly true since expanding the clump leaves it more susceptible to disruption through 
pressure and tidal stripping, a feature that would worsen the overmerging problem. 

7. DISCUSSION 

Using a series of simulations of a cold clump moving through a hot medium, we have explored 
the nature of drag in scenarios directly relevant to cosmological structure formation. The results 
are scale invariant and, subsequently, may be applied to astrophysical situations non-cosmological 
in nature. Indeed, any hydrodynamical simulation using the SPH technique in which drag on 
resolution-limited objects is significant will be affected by the results described herewith. 

An expression for a fiducial drag of the clump was determined analytically presuming com- 
pletely inelastic collisions. The expression derived, when supplemented by a coefficient of drag, 
is identical to previously derived analytic forms. The force goes as F^rag V^R^. An estimate 
was made for the cross-sectional radius of the clump which was found to be twice the averaged 
smoothing length, hij. This is an average of the SPH smoothing length common to the cold clump 
particles, hdumpi a-i^d that common to the hot gas particles, hhotj which is used to symmetrise the 
SPH force calculations. In the case of symmetrisation via kernel averaging, the radius of the clump 
should approximate 2hhot- 

The energy lost from the simulated clumps as they decelerate is deposited locally in the case 
of the subsonic clump. For a sonic and supersonic clump, SPH successfully models a bow shock 
which carries away the energy from the clump, depositing it over a larger volume. 

The drag was compared with the fiducial amount by equating the measured drag to the fiducial 
form and solving to determine an effective radius, Reff, for the clump. Over the entire range 

—3 / 

examined, the effective radius is fit quite well by the form, Reff o: M. . This implies that 

4/ 

Fdrag oc which is almost linear. 

In the regime M. >^/i-, the effective radii of the simulated clumps lie in the range Reff — hij 
to 2hij. This is consistent with the estimated value for the effective radius of 2hij, particularly if 
the coefficient of drag is between 0.5 and 1. However, there is no evidence of convergence to any 
particular value of Reff as A4 increases. 

For the range M. < ^U, the effective radius is larger than the estimate; i.e., Reff > 2hij. 
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Indeed, it is larger than even the maximum radius of direct hydrodynamic influence, 2hhot- This 
demonstrates that the drag on sub-sonic clumps exceeds any physically justifiable amount presum- 
ing inelastic collisions. 

It is notable that the transition point, A4 — ^A, is independent of the local velocity of sound. 
This demonstrates that the simulations are velocity-scale invariant. 

At all velocities, the drag force was found to scale with the cross-sectional area provided by 
the averaged smoothing length, hij. There arc two immediate implications of this: the method of 
force symmetrisation (which determines the method of averaging) affects the drag, and the effective 
cross section of the clump varies with the density of the medium through which it is travelling. 

In regards to the first implication, the method of force symmetrisation is selected during the 
coding of the SPH implementation. It can be one of the following four methods: arithmetically 
averaging the kernels, or averaging the smoothing radii of the interacting pair of particles using 
either arithmetic, harmonic, or kernel averaging. It was verified that harmonic averaging of the 
smoothing radii produces the smallest eff'ective cross-section whereas kernel averaging produces the 
largest. The eS'ect is substantial, with the drag produced by the kernel averaging 2.4 times larger 
than that of the harmonic averaging. The drag forces from geometric and arithmetic averaging of 
the smoothing lengths fall in between, with arithmetic averaging the larger of the two. 

Concerning the second implication of the cross-sectional radius scaling with hij, it was shown 
that there is a loss of proper scaling of the drag with density. The drag should scale with density 
as F oc p. This occurs because the average smoothing length, hij, varies with the density of the 
surrounding medium. However, higher densities lead to a smaller average smoothing length, which 
decreases the cross-sectional area of the clump, offsetting the increase in drag due to the increased 
density. This is particularly true for methods of force symmetrisation that use kernel averaging 
or arithmetic averaging of the smoothing lengths. For these, the drag on a clump is proportional 

1 / 2 / 

to p^^j instead of phot- For geometric averaging, this relation changes to Pf^^^, while for harmonic 
averaging, the effect is present over a limited range of hhot/hdump- 

This attenuation of the drag as the density increases will have consequences in simulations of 
the hierarchical formation of cosmological structure. A disproportionate amount of deceleration 
will occur in the outer halos of the structures. This will preferentially deposit energy in the outer 
halo as well as alter the spatial distribution of the deposited matter. This should directly affect 
the velocity and spatial bias of the baryonic material in the structures. 

More disturbing is the dramatically large drag felt by subsonic clumps. This was demonstrated 
to be a by-product of the method of force calculation used by SPH. For clumps travelling at 
Ai < '■'U, the method SPH uses to calculate the hydrodynamical forces, coupled with the compact 
nature of the cold clump, initiates a process that accretes hot particles onto a sphere about the cold 
clump. The hot gas particles are forced by the interaction with the cold clump to oscillate outside 
the radius provided by twice the averaged smoothing radius, triggering the artificial viscosity which 
boosts the velocity of the hot particles to match that of the cold clump. This process is inherent 
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to the SPH method since it arises from the force calculation. It leads to the excessive drag at 
low velocities because the particles accreted by the cold clump, or remoras, increase the effective 
cross-sectional area of the clump. 

The accretion of the remoras leads to an increase in drag by increasing the effective cross 
section. The flow of the particles in Fig. 5 indicates that the gas particles are not being repulsed 
solely by the clump, but also by the remoras in the shell about the clump. If this is the case, then 
the cross sectional radius should increase from 2hij to 2{hij + hhot)- For the systems presented 
here, this should increase the cross-sectional radius by 2.5 times. 

Comparing the effective radii calculated for the Mach Vs clump with those of the Mach 1, 
and Mach 2 clumps imply respective ratios of 2.3 ± 0.3 and 2.9 it 0.4 (see Table 2 for the radii). 
However, Fig. 3 shows that the change in i?e// does not occur as a step function, nor does it plateau 
at 2.5/ijj. This indicates that further processes are occurring. One possibility is 'crystallisation' 
of the chimp into the matrix of the hot gas particles as the velocity of the clump approaches the 
velocity dispersion of the clump (Lombardi, Sills, and Shapiro 1998). 

For virialised objects, the velocity dispersion of the component structures should be on the 
order of the speed of sound of the halo gas. This is approaching the velocity at which the accretion 
of remoras commences and we observe an inappropriate increase in the drag. This will enhance the 
merging rate of the structures. Not only will this modify the distribution of small-scale structures, 
but it, too, will affect the velocity and spatial bias of the baryons. 

Boosting the volume in which the cold clumps reside by providing a short-range anti-gravity 
force is shown to reduce the drag by a factor of two. This is not a viable solution to the problem, 
however, since the clumps are then not stable outside a hot halo. What may be required is a 
decoupling of the separate hot and cold phases, similar to that described in Pearce et al. (1999), 
with the inclusion of a solid-body force between the clump and the hot gas. 

We acknowledge NATO ORG 970081 which facilitated our interaction. We would like to thank 
Dr. R. Thacker for supplying the various implementations of SPH. This work was supported by 
NSERC Canada. 
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Slow cold clump Fast cold clump 



P/Pc 


1000 


1000 


T(K) 


10^ 


10^ 


R (kpc) 


50 


50 


N 


100 


100 


m (lO^M©) 


1.7 


1.7 


Vo (km/s) 


500 


1000 




Hot gas 


Very hot gas 


p/pc 


10 


10 


T(K) 


10^ 


10^ 


N 


13000 


13000 


m (10^ Mq ) 


1.7 


1.7 


Vs (km/s) 


500 


1500 


Rj (Mpc) 


6 


18 



Table 1: The characteristics of the cold clumps and the hot media used in the drag tests. Given are 
the overdensity, p/pc (^loo = 1)) the temperature, T, the radius of the cold clump, R, the number 
of particles in the medium, the mass resolution of the medium, m, the initial velocity of the cold 

clump, V^, the speed of sound in the hot medium, Vs, and the Jeans length for the hot medium, 
Rj. The simulation volume in all cases is (2Mpc)^. The 'fast cold clump' was used in the Mach 2 
runs in combination with the 'hot gas'. The Mach 1 runs used the 'slow cold clump' embedded in 
the 'hot gas'. The Mach 1/3 runs used the 'slow cold clump' in the 'very hot gas'. 
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Mach 1/3 



nhot 


i?e//(Mpc) 


Reff/ hclump 


Reff/hhot 


Reff/hij 


xl 


0.23 


6.0 


2.2 


3.3 


x2 


0.17 


4.4 


2.1 


2.8 


x5 


0.15 


4.0 


2.6 


3.1 


Mach 1 


nhot 


i?e//(Mpc) 


Ref f / hclump 


Reff/hhot 


Reff/hij 


xl 


0.13 


3.6 


1.3 


1.9 


x2 


0.11 


3.0 


1.4 


1.9 


x5 


0.07 


1.9 


1.2 


1.5 


Mach 2 


nhot 


Reff {Mpc) 


Reff 1 hclump 


Reff/hhot 


Reff/hij 


xl 


0.093 


2.5 


0.9 


1.3 


x2 


0.086 


2.3 


1.1 


1.5 


x5 


0.0(33 


1.7 


1.1 


1.3 



Table 2: Scaling of the effective radius, Reff, of the clump with mean smoothing length of the 
clump gas particles, hdump, the hot gas, hhot^ and the symmetrised smoothing length, hij, as used 
in the SPH calculations. Varied is the number density of the hot gas particles, nhot, but not the 
mass density. 
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Mach 1/3 



Phot 


^e//(Mpc) 


Reff 1 hclump 


Reff/hhot 


Reff/hij 


Xl 


0.25 


6.7 


2.5 


3.6 


x2 


0.22 


5.9 


2.8 


3.8 


x5 


0.19 


5.0 


3.2 


3.9 


Mach 1 


Phot 


i?e//(Mpc) 


R-eff/ hclump 


Reff/hhot 


Reff/hij 


Xl 


0.13 


3.6 


1.3 


1.9 


x2 


0.12 


3.1 


1.4 


2.0 


x5 


0.09 


2.4 


1.5 


1.9 


Mach 2 


Phot 


i?e//(Mpc) 


Reff/ hclump 


Reff/hhot 


Reff/hij 


xl 


0.097 


2.6 


1.0 


1.4 


x2 


0.089 


2.4 


1.1 


1.5 


x5 


0.068 


1.8 


1.2 


1.4 



Table 3: Scaling of the effective radius, Reff, of the clump with mean smoothing length of the 
clump gas particles, hdump, the hot gas, hhot^ and the symmetrised smoothing length, hij, as used 
in the SPH calculations. Varied is the physical density of the hot gas particles, phot- 
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0.5 1 1.5 2 2.5 3 3.5 4 




Fig. 1. — The kernel function for differing effective smoothing lengths. Plotted are the kernel 
functions for the symmetrisations of h via the arithmetic mean (large dots) and harmonic mean 
(small dots) when ^2 = 2^i. The kernel functions using just hi and /i2 are also given (dashed and 
dot-dashed, respectively). The effective kernel function for the average of these two kernels is also 
plotted (solid line). 
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Clump velocity (Mach) 



Fig. 2. — The mean drag felt on a clump moving with constant velocity. The data are given for 
both the very hot gas {Vs = 1500 km/s) and the hot gas {Vg = 500km/s). The lines indicate the 
theoretical maximum deceleration calculated in Sec. 5. 
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Clump velocity (Mach) 



Fig. 3. — The effective radius of the clump as it varies with Mach number. The data are given for 
both the very hot gas {Vg = 1500 km/s) and the hot gas {Vg = 500km/s). 



Fig. 4. — The temperature change in the hot gas due to the passage of the clump. The temperature 
changes for the gas is plotted for the Mach Va (top), Mach 1 (left), and Mach 2 clumps (right). 
White indicates a larger temperature increase, black indicates no change. Note that the temperature 
changes are small compared to the actual temperature of the gas (10^ K cf. 10^ K). The clump 
starts at a position one tenth of the box size from the left. The slices are lOhdump thick. 
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Fig. 5. — Flow of halo particles about clumps. The flows are plotted for the Mach Vs (top), Mach 
1 (left), and Mach 2 clumps (right). The inner circle is the radius 2hciump while the outer circle is 
the same for the hot gas, 2hhot- The effective smoothing length, hij, for this simulation is simply 
the arithmetic mean of hi and hj. Flow is from right to left. The slices are 5hciump thick. 
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Fig. 6. — The 2-point correlation function about the cold clump. Marked are the smoothing limits 
for the particles of the cold clump (dashed), the hot gas (dot-dashed), and the average (dotted) 
which is the value used in these tests for a hot particle interacting with the cold clump. For 
comparison, the 2-point correlation function for the hot gas away from the cold clump is given by 
the dashed curve. 
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X 

Fig. 7. — Path of a hot particle onto a cold, Mach Vs clump. The position of the particle at 
successive iterations is plotted relative to the cold clump. Flow is from the right (the clump is 
moving left). Again, the outer circle has a radius of 2hhot while the inner radius is 2hciump- 
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time (Ga) 



Fig. 8. — The forces on a remora particle that propel it in a direction perpendicular to the direction 
of travel. Plotted is the total force (solid line) as well as the contributions from the gas pressure 
(dashed line) and the artificial viscosity (dash-dotted line). The gas pressure provides the bulk of 
the instantaneous forcing and hence frequently overlaps the total force. Positive forces are directed 
toward the clump centre. At the bottom is plotted the gas density calculated by the SPH routine 
(dotted line). The axis to the right relates to this plot. 




Fig. 9. — The cumulative momentum gained by a remora particle due to SPH forces. The panel 
plots the forces gained in the direction of travel of the cold clump as it accumulated with time. 
Plotted is the total momentum boost (solid line) as well as the contributions of the SPH pressure 
term (dashed line) and the artificial viscosity (dash-dotted line). 



